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In recent publications, it has been shown that high-order harmonic generation can be manipulated 
by employing a time-delayed attosecond pulse train superposed to a strong, near-infrared laser field. 
It is an open question, however, which is the most adequate way to approximate the attosecond pulse 
train in a semi-analytic framework. Employing the Strong-Field Approximation and saddle-point 
methods, we make a detailed assessment of the spectra obtained by modeling the attosecond pulse 
train by either a monochromatic wave or a Dirac-Delta comb. These are the two extreme limits of a 
real train, which is composed by a finite set of harmonics. Specifically, in the monochromatic limit, 
we find the downhill and uphill sets of orbits reported in the literature, and analyze their influence 
on the high-harmonic spectra. We show that, in principle, the downhill trajectories lead to stronger 
harmonics, and pronounced enhancements in the low-plateau region. These features are analyzed 
in terms of quantum interference effects between pairs of quantum orbits, and compared to those 
obtained in the Dirac-Delta limit. 



I. INTRODUCTION 

One attosecond (10 _18 s) is roughly the time it takes for 
light to travel through atomic distances. This fact makes 
high-frequency pulses of attosecond duration a very pow- 
erful tool for resolving or even controlling dynamic pro- 
cesses occurring at the atomic scale Indeed, in the 
past few years, such pulses have caused a breakthrough in 
metrology, with applications as diverse as resolving the 
motion of bound electrons 0, exciting inner-shell elec- 
trons 3|, inducing vibrational wavepackets in molecules 
0, or controlling electron emission 

Attosecond pulses have been predicted theoretically in 
the mid-1990s and have become experimentally feasi- 
ble a few years later 0, & 13 • They are a direct conse- 
quence of the fact that high-order harmonics, generated 
by the interaction of a strong laser field (intensities of 
the order of 10 14 W/cm 2 ) with a gaseous sample, are al- 
most phase locked [6| . Hence, by superposing sets of high 
harmonics, it is possible to produce attosecond pulses. 
Specifically, there exist two main scenarios: either one 
obtains an attosecond-pulse train, by employing long in- 
put pulses or isolated attosecond pulses, from few- 
cycle driving fields |2|, la [l(J . Both situations have been 
widely exploited in the literature, and sometimes the at- 
tosecond pulses appear in combination with additional 
driving fields. 

For instance, recently, an attosecond-pulse train su- 
perposed to a strong, infrared, linearly polarized field 
has been em ploy ed to control high-harmonic generation 
(HHG) 0] and above-threshold ionization (ATI) 
|13| . Such a train exhibited a time-delay td with respect 
to the infra-red field. By varying this delay, one could 
influence several features, such as the intensities, reso- 
lutions, and maximal energies in the high-harmonic or 
photoelectron spectra. This scheme has been proposed 
theoretically 0,0] and realized experimentally [l2 . ll3f . 



and the key idea behind it is to provide an additional 
pathway, which can be controlled, for an electronic wave 
packet to be released in the continuum. 

This can be understood in view of the physical 
mechanisms behind both phenomena. At a time 
t' , an electron is freed by tunneling or multiphoton 
ionization. Subsequently, it propagates in the continuum, 
gaining energy from the field, and it is driven back to- 
wards its parent ion, to which it returns at a time t. If 
the electron recombines or rescatters, there is either emis- 
sion of high-frequency radiation (i.e., HHG), or of high- 
energy photoelectrons (i.e., high-order ATI) 0j. For 
low-order ATI peaks, the electron reaches the detector 
without rescattering. An attosecond pulse train allows 
the electron to reach the continuum by absorbing high- 
energy photons and being able to overcome the ionization 
potential. Additionally, since they are of very short du- 
ration, the attosecond pulses provide a tool to control the 
instant at which the electron is being ejected. This has 
direct consequences in the momentum it acquires from 
the field when it is released 0], and also a strong influ- 
ence on the kinetic energy the electron has upon return 
0,0]. Therefore, it also affects the spectra. 

The results in 0,0] and 0] have been obtained from 
the numerical solution of the time-dependent Schrodinger 
equation, and the attosecond-pulse train has been mod- 
eled by the superposition of five harmonics (fin = 11a; 
to fiig = 19uj). These are realistic assumptions, and well 
within the available experimental data. However, in or- 
der to interpret the results obtained, the attosecond-pulse 
train has been approximated by a monochromatic wave 
of frequency firs = 15w, which is the central frequency of 
the group of harmonics in question, and a classical model 
for an electron ejected by such a wave and propagating in 
a strong, low-frequency field was constructed. These sim- 
plifications have enormously facilitated the understand- 
ing of the ionization mechanisms, and, specifically, have 
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shown how to manipulate ionization by varying the time 
delay t d . 

In a previous publication, we have also investigated 
the ATI and HHG spectra from an atom irradiated by 
a time-delayed attoseco nd-p ulse train superposed to a 
near-infrared laser field 15]. Specifically, we have con- 
sidered the transition amplitudes for both phenomena 
within the Strong-Field Approximation, and took the at- 
tosecond pulse train to be a sum of Dirac-Delta functions 
in the time domain. These assumptions have allowed an 
almost entirely analytical treatment, which has been vital 
for a clear understanding of the problem. Indeed, we were 
able to investigate the influence of the attosecond pulses 
on the spectra in far more detail than in [TJ fl2l Il3| , and 
to interpret the results obtained in terms of quantum in- 
terference effects. 

To a very large extent, our results, as well as their 
physical interpretation, agree with those in [TTL IT2I Il3| . 
In fact, we have obtained, for very specific time delays, 
enhancements of more than one order of magnitude in 
the harmonic spectra similar to those reported in |ll| . 
Furthermore, as the time delay increases, the enhance- 
ments move gradually from the high-order harmonics to- 
wards lower frequencies, until they only affect the low- 
plateau region. Finally, for ATI, we have observed an 
identical behavior of the yield with the time delay as in 
[l3j] . namely that it extends towards higher energies if 
ujtd = nir, and that it decays most quickly, as the fre- 
quency increases, for u>td = (2n + 1)tt/2. 

There exist, however, a few discrepancies, especially 
for HHG, between our results and those in [ll|. First, 
the above-mentioned features occur for different phases 
4> = utd- Second, the explanations for such features are 
slightly different. In we related the enhancements 
in the spectra to a particular pair of orbits, for which 
the excursion times At — t — t' of the electron in the 
continuum were very short. Therefore, the spreading of 
the electronic wave packet in the continuum would be 
very reduced. Hence, the overlap of this wave packet, 
upon return, with the bound state it left, would be con- 
siderable, leading to very prominent harmonics. We have 
shown that the maximal kinetic energy the electron may 
have upon return, for this pair of orbits, was very much 
dependent on the delay between the infra-red field and 
attosecond-pulse train, going from 1.8U p , where U p is the 
ponderomotive energy, to vanishingly small values. In 
the former case, this leads to strong harmonics through- 
out, whereas the latter case causes enhancements only in 
low-energy regions of the spectra. 

In Ref. [ll|], however, the modifications in the spectra 
have been attributed to the selection of specific electron 
trajectories, which are present in case it is released by 
a high-frequency wave, and, subsequently, propagates in 
a strong infra-red field. In comparison to the case in 
which the electron is released only by the infrared field, 
there is a splitting in the electron orbits, into the so- 
called "downhill" and "uphill" trajectories, with respect 
to the effective potential barrier V e g = V(r) — r.E(i') 



at the electron start time t', where E(i') and V(r) de- 
note the electric field and the atomic binding potential, 
respectively. The former and the latter case relate to 
the situation for which the electron velocity, at t', has 
the opposite or the same direction of the field, respec- 
tively. If the time delays are appropriately chosen, one 
may enhance ionization, and thus the HHG spectra, for 
a particular set of orbits, or for none. 

Probably, such discrepancies are rooted on the fact 
that the attosecond pulse train has been modeled in dif- 
ferent ways. While in I 111 it has been approximated by a 
monochromatic wave, in |l5j we considered a Dirac-Delta 
comb. In this context, one should keep in mind that a re- 
alistic attosecond-pulse train possesses a broad frequency 
range. Therefore, a monochromatic wave is as much as 
an approximation as a Dirac-Delta comb. The main dif- 
ference lies in the number of harmonics composing the 
train: the former should work better if this number is 
small (such as in |ll|), and the latter if it is large. 

In this proceeding, we present a detailed discussion of 
the two limits which one may adopt, when modeling the 
attosecond pulse train. We put a particular emphasis 
on the consequences of taking the attosecond pulse train 
to be a monochromatic wave, within the Strong-Field 
Approximation, for high-order harmonic generation. In 
order to make an assessment of the similarities and dif- 
ferences from the Dirac Delta comb employed in our pre- 
vious work ^5| , we will also briefly recall such results. 

This manuscript is organized as follows: In Sec. [nj 
we present the explicit expressions for the transition am- 
plitudes, if the attosecond pulse train is taken to be a 
monochromatic wave (Sec. Ill All or a Dirac Delta comb 
(Sec. IIIBjl . In particular, we discuss how the saddle- 
point equations are modified for each case, and relate 
their solutions to the classical orbits of an electron in a 
field. Subsequently, in Sec. IIIII we compute such orbits 
explicitly, and employ them to calculate high-harmonic 
spectra. Finally, in Sec. IIVI we conclude the paper, re- 
lating both limits to the results in [TT| . 

II. TRANSITION AMPLITUDES 

The general expression for the HHG transition ampli- 
tude, in the strong-field approximation |18|. is given by 

/■OO ft /> 

M Q = i dt dt' I d 3 kd* z (k + A{t))d z (k + A(i')) 

J — ocJ — OO J 

E{t')exp[iS(t,t',Q,k)}, (1) 
with the action 

1 r* 

s(t,t\n,k) = -- / [k+A{T)] 2 dT-i p {t~t')+nt, (2) 

2 Jr 

where A(t), J7,k and I p denote the vector potential, the 
harmonic frequency, the intermediate momentum of the 
freed electron and the atomic ionization potential, respec- 
tively 0] . Eq. (JTJ describes a physical process in which 
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an electron is freed at a time t', propagates in the con- 
tinuum with momentum k from t' to t, and, at this time, 
recombines with its parent ion, generating a harmonic of 
frequency f2. The prefactors 

d(k + A(t)) = (i Jo \r.e x \k + A(t)) (3) 

and 

d*(k + A(t')) = <k + A(f)| v.e x |Vo) , (4) 

where Vo) and E{t') = —dA(t')/dt' denote the electronic 
bound state and the electric field at the time the electron 
is released, respectively, contain all the influence of the 
atomic binding potential, which is implicit in the bound- 
state wave functions. In general, the vector potential is 
given by A(t) = A;(t) + A/j(i), and the total electric field 
by E(i) = E;(t) + E/^t), where the indices I and h refer 
to the laser field and the attosecond pulses, respectively. 
In this paper, we work within the so-called "broad gaus- 
sian limit" , which implies that the electron is bound, and 
returns to, a bound state localized at the origin of the co- 
ordinate system [l^. This yields constant prefactors © 
and (@J). 

For low enough frequencies and high enough laser in- 
tensities, Eq. Q can be solved to a good approximation 
by the steepest descent method. Thus, we must deter- 
mine k, t' and t so that S(t, t', Q, k) is stationary, i.e., so 
that its partial derivatives with respect to these param- 
eters vanish. Apart from considerably simplifying the 
computations involved, this approximation has the ad- 
vantage of providing a clear space-time picture for the 
physical process in question, and allowing a detailed as- 
sessment of quantum-interference effects. Furthermore, 
as it will be shown below, the saddle-point equations can 
be related to those describing the classical orbits of an 
electron 

If the attosecond pulses are absent, the saddle-point 
equations read 

[k + A l {t')f = -2I p (5) 



2(fi-I p ) = [k + A,(t)] 2 , (6) 

and 

Jdr[k + Ai( T )]=0. (7) 

Eq. JSJ and J§} give the energy conservation at the 
start and recombination times, respectively, while Eq. 
(|7|) yields the intermediate electron momentum fulfill- 
ing the condition for the electron to return. The first 
equation expresses tunneling ionization at t', and has no 
real solution. Physically, this means that tunneling has 
no classical counterpart. In the limit of I p — > 0, this 
equation describes a classical electron being released with 
vanishing drift velocity. Eq. © illustrates a process in 



which the kinetic energy of the electron upon return is 
converted into a photon of frequency Q. Near the cutoff, 
whose energy position corresponds to the maximal value 
this quantity may take, such a process has no classical 
counterpart, and the transition probability is expected 
to decay exponentially. 

We will now include the attosecond-pulse train. The 
strong, low-frequency field will be approximated by a lin- 
early polarized monochromatic wave of amplitude Eq, 
i.e., 

E l (t) = E sm(ujt)e x , (8) 
while the attosecond-pulse train will be taken as 

E 2*i+l 

Eh(i) = ^r V wx[qu{t-t d )]e x , (9) 
cr(t) 

v ' <;=2fco + l 

of amplitude E^. Eq. 10 is a sum over several odd high- 
order harmonics, which are phase-locked and exhibit a 
time delay td with respect to the background field. The 
indices 2ko + 1 and 2k\ + 1 yield the minimal and the 
maximal harmonic order, respectively. The function a(t) 
denotes the train temporal envelope, which will be taken 
as (7 = const. Physically, this corresponds to an infinitely 
long attosecond-pulse train. 

In our model, we assume that the attosecond pulses 
will mainly influence the electron ejection in the con- 
tinuum but not its subsequent propagation. This latter 
step will be governed by the low-frequency field. For this 
reason, unless stated otherwise, we will use the approx- 
imation E(t') ~ Eh(t') in the prefactor of Eq. QJ, but 
not in the action. This implies that the vector potential 
is approximated by A(t) ~ Ai(t). 



A. The monochromatic limit 

For a finite number of harmonics, this yields highly 
oscillating prefactors, which are no longer slowly varying. 
Hence, they must be incorporated in the action. The 
HHG transition amplitude then reads 

2fci+l 

M » = -|| E M t q \ do) 

< ? =2feo+l 



with 

/oo ft r 
dt I dt' d 3 kd*(k + Ai(t))d(k + A/0')) 
-oo J — oo J 

^(t'-u) _ e -n"(t'-U)) exp [iS(t, t', n, k)](ll) 
In Eq. (|10|l . we will define a modified action, 

S = S±qu(t' -t d ). (12) 
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This leads to changes in the saddle-point equation JSJ, 
which now reads 



[k + Mt')} 2 = 2(=f«w - /„)• 



(13) 



Thereby, the solution in — quj does not make physically 
sense, and thus will not be taken into account, whereas 
that in +quj describes an ionization process in which an 
electron is ejected by the absorption of a high-frequency 
photon. Hence, the attosecond-pulse train, in which a 
superposition of such processes occurs, provides an alter- 
native pathway for the electron to reach the continuum. 
This is in agreement with the discussions in ^jj . We will 
consider here the limit for which the attosecond pulse 
train is described by a monochromatic wave, so that the 
sum in 1)10(1 is dropped. 

If Eq. (|13fl is combined with JBJ and J7J) , one obtains 
the expressions 



ei arccosa 



(14) 



and 



£2vl — a 2 — smut = (15) 
= u>(e± arccosa — t)(cosut — 72), 

with a = cosut + 71 — 72, and = ±1 (i = 1,2), for 
the electron ionization and return times, respectively. 
Hence, the problem of finding the pairs of times (t, t 1 ) 
has been reduced to solving the transcendental equation 
(JT5J. In the above-stated equations, one may identify the 
parameters 



7i 



±4 



IQ-L 



2U 7 



- and 72 = ±i 



2U„ 



(16) 



In particular 72 is very important for determining the 
initial conditions for the electron when it is being ejected 
in the continuum. 

One should note that there exist sign ambiguities in the 
above-stated equations. Furthermore, additional com- 
plications may be caused by the fact that arccosa = 
arccos(2rt7T — a). We will discuss such problems explicitly 
in Sec. IIII Al in which they are overcome by means of 
physical arguments. 



B. The Dirac-Delta limit 

In the limit k\ — ko — > 00, which is the main assump- 
tion in our former work |l5j |. Eq. ® reads 



E h (t) 



J2(-l) n S(u ) (t-t d )-n7r)e x , (17) 



n=0 



This yields the transition amplitude 



(D) _ 



^Hi- l ) n J dt J d 3 fcexp[z5(t,4,^,k)] 



In this case, the ionization time is being fixed at t' n — 
td + nir/u>. The above-stated equation means, physically, 
that the electron is no longer reaching the continuum 
through tunneling ionization, but is being released by the 
attosecond pulses. However, in contrast to the approxi- 
mation in the previous section, all information about the 
velocity with which the electron is leaving is lost. In- 
deed, the electron is reaching the continuum with any of 
the energies Nu — I p , since all harmonics composing the 
train are equivalent. 

Under these assumptions, the transition amplitude is 
further simplified so that only the integrals in the in- 
termediate electron momentum k and the return time t 
must be solved. In this case, the saddle-point equations 
iJSJ and J7J can be combined in order to obtain the tran- 
scendental equation 



smut — (—1)™ smutd 

[u(t — td) — nir] (cosut — 71) 



(19) 



^(k + A,(t))d,(k + A,(0). 



(18) 



where 71 is defined in l|16fl . Similarly to the 
monochromatic-limit case, one may identify ambiguities 
in Eq. I|19l) . It turned out, however, that they are fewer 
and far easier to overcome in the limit discussed in this 
Section. 



III. RESULTS 
A. Start and return times 

We will now assess the influence of the additional 
attosecond-pulse train on the electron orbits. As a start- 
ing point, we will approximate it by a high-frequency 
monochromatic wave. In general, we will refer to a pair 
of orbits employing the indices which increase with 

the electron excursion times At = t — t' in the continuum. 

We will initially concentrate on the shortest pair of or- 
bits, i.e., (1,2), which provide the most prominent con- 
tributions to the high- harmonic spectra. In Fig. 1, we 
display the real parts of the ionization and return times 
for such orbits, obtained from the saddle point equations 
1 31) , I© and J7J) , as compared to the situation for which 
the attosecond pulses are absent (i.e., for q = in IfT^i ). 
In all cases, start and return times occur in pairs, which 
coalesce at well-defined energies. Such energies corre- 
spond to the maximal kinetic energy with which a clas- 
sical electron may return to its parent ion. For higher 
energies, the transition amplitudes are exponentially de- 
caying. Hence, potentially, such energies lead to cutoffs 
in the spectra. For a monochromatic field in the ab- 
sence of the attosecond-pulse train, the cutoff occurs at 
3.17Up + I p . 

In the figure, we can identify two main distinct regimes. 
If the frequency Q q = qui is lower than the ionization po- 
tential I p , qualitatively, there is the same behavior for 
the start times t' or return times t. Indeed, the only dif- 
ference is the cutoff energy, which is slightly lower. An 
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FIG. 2: (Color online) Imaginary parts of electron start times 
for the two shortest orbits of an electron subject to the same 
field and atomic parameters as in the previous figure. In the 
figure, we only show the harmonic orders for which f2 > I p . 
The dashed and the solid lines in the figure refer to the orbit 
1 and 2, respectively. The numbers on the upper edge of the 
figure denote the kinetic energy of the electron upon return, 
in units of the ponderomotive energy U p . 



FIG. 1: (Color Online) Real parts of return and start times 
(panel (a) and (b), respectively), as functions of the high- 
harmonic order, for the two shortest orbits of an electron sub- 
ject to a monochromatic field of intensity I = 4 x I0 14 W/cm 2 , 
and frequency to — 0.057 a.u., superposed to a high-frequency, 
monochromatic wave of Ih = I/W, approximated by a 
monochromatic wave of frequency Q q = qui. For compar- 
ison, we also provide such times for the case in which the 
attosecond train is absent (q — in Ea. l|I3^ 'l. The atomic 
system is taken to be argon (I p = 0.58 a.u.) and we only 
show the harmonic orders for which Q > I v . The dashed and 
the solid lines in the figure refer to the orbit 1 and 2, respec- 
tively. The numbers on the upper edge of the figure denote 
the kinetic energy of the electron upon return, in units of the 
ponderomotive energy U p . 



inspection of Eqs. (|14|1 and l|15Jl supports this inter- 
pretation: if quj — - I p < 0, the parameter 72 is purely 
imaginary, so that the solutions of (|15|) occur in conju- 
gate pairs. The element of such a pair which leads to 
Im[i'] < is discarded as unphysical, whereas the other 
element is directly associated with the tunneling process. 
In this context, the quantity Im[£'] provides information 
about how probable the tunneling process is, as it will be 
discussed below. 

If on the other hand, qtu — I p > (for instance, for 
q = 15), the attosecond pulses have given enough en- 
ergy for the electron to overcome the potential barrier. 
In this case, it is reaching the continuum with veloc- 
ity v = ±y/2(quj — Ip). In comparison to a situation for 
which the electron is being released with vanishing ve- 



locity, which is the classical limit of Eq. (0, there is a 
splitting in the ionization times t 1 for each set of orbits. 
This is precisely the effect reported in Ref . ^l| , and can 
be readily understood from Eq. I|15|) . in which now the 
parameter 72 is real. Hence, there exist now two solutions 
of (|15|) which make physically sense. The solution with 
72 < yields the downhill trajectories in 01, whereas 
that with 72 > corresponds to the uphill trajectories. 

Interestingly, however, for both types of trajectories, 
the electron returns almost at the same time t. In fact, 
there are only noticeable differences near the ionization 
threshold Vl = 10w, which, at most, will influence the 
low-plateau region. Hence, the electron excursion times 
in the continuum are considerably shorter for 72 < 0. 
Obviously, the larger the frequency fl q = qtu, the more 
pronounced the splitting is. In comparison with 
however, there is a much larger overlap of the start times 
for both types of trajectories near the ionization thresh- 
old. 

From a more technical viewpoint, it is worth mention- 
ing that, in the computation of the orbits (1, 2), we have 
taken ei = £2 = +1 and 71 < in general. As an ex- 
ception, for q = 15, 72 < 0, and i = 1, for harmonic 
frequencies f2/w < 24.4 one must employ an analytical 
continuation in Eq. (|15fl . with 62 = —1, £i = +1 and 
27T — arccosa instead of arccosa, in order to guarantee 
continuous solutions. 

The above-stated physical interpretation is confirmed 
by Fig. 2, in which the imaginary parts of the start 
times are depicted. Such imaginary parts are in a sense 
a measurement of a quantum-mechanical process having 
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FIG. 3: (Color online) Real parts of start and return times 
for the pairs of orbits (3,4) [panels (a) and (b)], and (5,6) 
[panels (c) and (d)], as functions of the high-harmonic order, 
for the same atomic and field parameters as in the previous 
figures. The orbits (3, 4) have been obtained by taking ei = 
£2 = +1 and 71 > in Eq. I15H . while (5,6) have been 
found employing e± — e-2 — +1 and 71 < 0. The numbers on 
the upper edge of the figure denote the kinetic energy of the 
electron upon return, in units of the ponderomotive energy 



a classical counterpart. Indeed, the smaller |Im[i']| is, the 
larger the probability that the process in question takes 
place will be. In the figure, we observe that, for the 
monochromatic and q = 9 cases, |Im[t']| 7^ 0. This is ex- 
pected, since tunneling ionization has no classical coun- 
terpart. Furthermore, |Im[t']| decreases if the attosecond 
pulses are present, due to an effectively narrower barrier. 
In case the electron is able to overcome the atomic bind- 
ing potential (e.g., for q = 15), for energies lower than 
the cutoff, |Im[t']| is vanishingly small, since now a clas- 
sical counterpart does exist. Clearly, beyond the cutoff 
energy, the imaginary parts of such times increase. 

The splitting in the electron start times is not exclusive 
of (1,2), but, in fact, occurs for all pairs of orbits. As 
concrete examples, in Fig. 3 we illustrate the real parts 
of the start and return times for the orbits (3, 4) and 
(5, 6). Such orbits correspond to electron excursion times 
At ~ 1.5T and At ~ 2T, respectively, where T = 2-k/uo 
denotes a cycle of the driving field. Clearly, the start 
times t' split for qui — I p > 0, due to the electron non- 
vanishing velocity upon ejection. Once more, each pair 
of start times leads to a shorter and a longer set of orbits, 
which correspond to 72 < or 72 > 0, respectively. The 
corresponding electron start times do not overlap in any 
energy region. This is due to the fact that they are much 
more localized than for the orbits (1,2). 

For comparison, in Fig. 3, we present the solutions 
of Eq. Q19f) . obtained when the attosecond-pulses are ap- 



FIG. 4: (Color online) Real part of the recombination times 
as a function of the high-harmonic order, for neon (I p — 0.79 
a.u.) interacting with a monochromatic field of intensity 
I = 5 x 10 14 W/cm 2 , and frequency u = 0.057 a.u. and an 
attosecond-pulse train of intensity Jh = 10 13 W/cm 2 , com- 
posed of an infinite number of harmonics (Eq. 117^ 1. Parts 
(a), (b), (c) and (d) correspond to delays = 0, 4> = 0.25-7T, 
4> — 0.5-7T and <f> = 0.75-7T between the attosecond-pulse train 
and the low-frequency driving wave, respectively. In the fig- 
ure, we only show the harmonic orders for which fl > I p , and 
the pairs of orbits are indicated by the natural numbers 
which increase according to the electron excursion time in the 
continuum. The ionization times are indicated by the blue 
dashed lines, and the numbers on the upper edge of the fig- 
ure give the approximate kinetic energy of the electron upon 
return, in units of the ponderomotive energy U p . 

proximated by a sum of Dirac-Delta functions, for several 
time delays. The start and return times, as well as the 
cutoff energies, are different from those in the previous 
figures. Indeed, we no longer observe a splitting in the 
sets of ionization and recombination times, as compared 
to the situation for which the attosecond pulses are ab- 
sent. This is in perfect agreement with the fact that now 
there is no longer a well-defined velocity with which the 
electron is being ejected. In contrast, however, there is 
now a single ionization time per half cycle, determined 
by Eq. I|17|) . The return times, however, still occur in 
pairs that coalesce at the cutoffs. For each pair of or- 
bits, the cutoff energies, as well as the excursion times, 
considerably changes with the time delay td. 



B. Harmonic spectra 

We will now analyze the spectra computed under the 
assumptions in Sec. Ill Al i.e., that the attosecond- 
pulse train can be approximated by a high-frequency, 
monochromatic wave. For comparison, we will also 
present spectra computed considering that the APT is 
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a sequence of Dirac-Delta functions. In this latter case, 
however, we will keep the discussion as brief as possible, 
as a detailed analysis is already presented in . In both 
cases, we have employed a uniform saddle-point approx- 
imation in order to compute the transition amplitudes. 
This approximation is discussed in detail in 

If one concentrates on the physical picture of down- 
hill and uphill trajectories, a very important question is 
the influence of each type of orbit in the high-harmonic 
spectra, and how such trajectories can be selected. In 
such a selection has been performed by employing 
an appropriate time-delay td between the strong, infra- 
red field and the attosecond-pulse train. In our frame- 
work, however, since we are approximating the train by a 
monochromatic wave, all information about this param- 
eter is lost. However, we can mimic the effect achieved 
in [Tll | by leaving out, or including, the corresponding 
orbits in our computations. 

In Fig. 5. (a), we present the spectra computed for a 
train with the frequency Q15 = 15w employing all orbits 
up to (5, 6). In this case, quo~I p > 0, so that the electron 
is being ejected with non-vanishing drift velocity and the 
spectra contains contributions from both the uphill and 
the downhill trajectories. We have restricted the electron 
start times to the first half cycle of the laser field, so that 
there are no sharp harmonics in the figure. One is able to 
view, however, the structures caused by the interference 
between the different pairs of trajectories very clearly. 
Prominent features are also a cutoff near harmonic order 
N = 57. This is in agreement with Fig. 1, which shows 
that the pairs of orbits (1,2) coalesce around this energy, 
and also with Fig. 2, for which there is a sharp increase 
in |Im[i']| at this harmonic order. Furthermore, there is 
a very pronounced hump in the low-plateau region. 

More detail is provided in Fig. 5.(b), in which the indi- 
vidual contributions from the downhill (72 < 0) and the 
uphill (72 > 0) trajectories are displayed. In general, the 
former contributions are at least half an order of magni- 
tude larger than the latter. This is expected, since, for 
72 < 0, the electron excursion times in the continuum 
are shorter for all sets of orbits. This means that there 
is less spreading for the electronic wave packet, and con- 
sequently a larger overlap, upon return, between such a 
wave packet and the atomic bound state with which it 
recombines. Obviously, this leads to a larger harmonic 
yield. Furthermore, in both cases there is no difference 
in the cutoff energy. This is due to the fact that all tra- 
jectories split symmetrically, with respect to the purely 
monochromatic case. 

The above-stated features can be understood in detail 
analyzing the contributions from each pair of orbits. Such 
contributions are displayed in Fig. 6. For both cases, the 
main features in the spectra, such as their overall inten- 
sities, the plateau shapes and cutoff energies, are deter- 
mined by the orbits (1, 2). In fact, the contributions from 
such orbits are at least one order of magnitude larger 
than those from (3, 4) or (5, 6). This is in particular true 
for the case 72 < 0, due to the shorter excursion times. 




Harmonic Order N 



FIG. 5: Harmonic spectra computed with the orbits (1,2), 
(3,4) and (5,6), for an attosecond-pulse train of frequency 
0,15 = 15a; superposed to a strong, near-infrared field of fre- 
quency u — 0.057 a.u. The field and atomic parameters are 
the same as in Fig. 1. Panel (a) exhibits the contributions 
from all sets of orbits, while panel (b) displays the individ- 
ual contributions from the orbits with 72 > or 72 < in 
Eq. 11511 . The numbers on the upper edge of the figure give 
the kinetic energy of the electron upon return, in units of the 
ponderomotive energy U p . 



Apart from that, they coalesce at a much larger energy 
than the other sets. This means that, classically, the 
maximal energy for which the electron returns to its par- 
ent ion along (1,2), and consequently the cutoff, is larger 
than for the longer orbits. In the figure, we also notice 
that the hump in the low-plateau region, for 72 < 0, 
can be traced back to the fact that the electron excur- 
sion time in the continuum, for the set (1, 2) is extremely 
short in this case (see, e.g., Fig. 1). The other pairs of 
orbits mainly contribute to the substructure in the spec- 
tra. Another noteworthy feature is that, regardless of the 
sign of 72 , the cutoff energies are the same for all sets of 
orbits, in agreement with the previous discussion. 

Another scenario, presented in Fig. 7 for comparison, 
is observed if the attosecond-pulse train is taken to be 
a Dirac-Delta comb. In this case, the dominant set of 
orbits, and, consequently, the shape of the spectra, are 
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FIG. 6: Individual contributions from the shortest three pairs 
of orbits to the high-harmonic spectra, for the same field and 
atomic parameters as in the previous figure. Panel (a) and 
(b) depicts the contributions from 72 > or 72 < in Eq. 
1151 . respectively. For comparison the overall contributions 
are displayed as the solid black lines in the figure. The num- 
bers on the upper edge of the figure give the kinetic energy 
of the electron upon return, in units of the ponderomotive 
energy U p . 



highly dependent on the time delay td- For instance, if 
the time delay is vanishing [Fig. 7. (a)], the spectrum is 
dominated by the orbits (1,2). This is due to the fact 
that they coalesce at the energy I p + 1.8U p , which is at 
a relatively high harmonic order. Furthermore, since the 
electron excursion times are very short, the harmonics are 
quite prominent. For <j> = utd = 0.257T [Fig. 7.(b)], there 
is now a double plateau, with two consecutive intensity 



drops at I p + 0.93U P and I p + 2.11U p . These are the 
energies for which the orbits (1,2) and (3,4) coalesce, 
respectively. Specifically in the upper half of the plateau, 
there was a drop in one order of magnitude, as compared 
to Fig. 7. (a). This is due to the larger excursion times, 
and, consequently, wave-packet spreading, for (3,4). If 
the delay increases further [Fig. 7.(c)], the contributions 
from (1, 2) only cause a shoulder in the spectrum, due to 
the fact that they now coalesce at I p + Q.2QU P , and the 
plateau, extending until I p + 2.5U P , is determined by the 
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FIG. 7: (Color online) Contributions of isolated pairs of 
orbits to the high-harmonic spectra for the same atomic 
and external-field parameters in Fig. 4, and delays cj> = 
0, (j> — 0.25-7T, 4> — 0.57T and (f> = 0.757T between the 
attosecond-pulse train and the low-frequency driving wave 
(Figs- Eta), [71(b), |2(c) and|2(d), respectively). In the fig- 
ure, we only show the harmonic orders for which Q. > I p . The 
contributions from specific pairs of orbits are depicted in the 
same colors as those employed in Fig. 4, for depicting each 
corresponding pair. For comparison, the spectra computed 
with the ten shortest orbits and with orbit 3 to 10 are given 
as the thick lines in the figure. The numbers on the upper 
edges of the figure give the approximate kinetic energy of the 
electron upon return, in units of the ponderomotive energy 
U v . 



orbits (3,4). The latter orbits dominate the spectra even 
more for <fi = utd = 0.757T [Fig. 7.(d)], since the cutoff 
energy for (1,2) is now vanishingly small. This effect 
is very similar to that observed in |ll|. even though it 
occurs for completely different delays. 



IV. CONCLUSIONS 

The outcome of this work supports the results in [Tl| in 
several ways. In general, it also shows that, by changing 
the initial conditions with which an electron is ejected 
in the continuum, one may manipulate several features 
in the spectra, such as the high-harmonic intensities, or 
cutoff. Furthermore, by considering the electron to be 
released in the continuum by a monochromatic wave, we 
have obtained a splitting of the electron start times very 
similar to that reported in |ll| , within the context of the 
Strong-Field Approximation, in the so-called uphill and 
downhill trajectories. 

We have also observed enhancements in the low- 
plateau region, and higher overall intensities for the con- 
tributions from the downhill trajectories. This indicates 
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a good perspective for harmonic control, by selecting ei- 
ther the downhill or uphill trajectories with, for instance, 
an adequate delay between the infrared field and the 
attosecond-pulse train. 

We have also obtained very detailed information about 
the contributions of specific pairs of orbits to the spectra. 
In the monochromatic limit, the shortest pair of trajecto- 
ries, denoted (1,2) is completely dominant, determining 
the shape of the plateau, the cutoff energy and the over- 
all intensity in the spectra. This is due to the fact that 
the maximal classical energy for this pair, and thus the 
cutoff, is always at a higher harmonic order than for the 
longer pairs of trajectories. 

In contrast, in case one considers a sum of Dirac-Delta 
functions, the dominance of this pair will depend very 
much on the time delay between the infra-red field and 
the attosecond pulses. This is due to the fact that the 
maximal energy an electron returning along these orbits 
is not always at a higher harmonic frequency than for 
the longer pairs. For instance, if for ujtd = 0, the or- 
bits (1,2) coalesce at the end of the plateau, and thus 
overshadow the remaining contributions, for tutd = 0.257T 
and Lotd = 0.57T, this energy lies at the low-energy part of 
the plateau, or near the ionization potential, respectively. 
Hence, the orbits (3,4) will also play an important role 
in the spectra. 

A feature in common for both the monochromatic 
and the Dirac-Delta limits is that the excursion time 
At — t — t' of the electron in the continuum is a very im- 
portant parameter, and considerably influences the shape 
of the plateau and harmonic intensities. This is not sur- 
prising, since, now, for all sets of orbits, the electron is 
reaching the continuum with a roughly equal, and large, 
probability. Hence, the differences in the yield will be 
mainly determined by the overlap between the returning 
electronic wave packet and the bound state with which it 
recombines. The shorter the orbit along which the elec- 
tron is returning is, the larger this overlap will be. In 
contrast, if the electron is released by tunneling ioniza- 



tion, the probability with which it is ejected will depend 
much more critically on the instantaneous potential bar- 
rier, and also influence the yield. 

Moreover, it is remarkable that both asymptotic limits 
exhibit features in common with those obtained in p"H 
within a more realistic model for the attosecond pulses. 
The monochromatic limit, for instance, allows one to 
identify the downhill and uphill trajectories, and to an- 
alyze their consequences in the spectra. Approximating 
the attosecond pulse train by a sum of Dirac Delta func- 
tions, on the other hand, sheds light on several features 
reported in such as a double plateau, or changes of 
more than one order of magnitude in the harmonic spec- 
tra. However, only in Delta-Dirac case there is a decrease 
in the energy of the enhanced harmonics, with increasing 
time delays td- 

There exist, however, aspects of ^lj that we have not 
addressed. We have not, for instance, investigated the 
differences in the high-harmonic resolution, due to a par- 
ticular trajectory choice. Indeed, in order to do so, it 
would have been necessary to extend the start times to 
several cycles of the laser field. Instead, we have re- 
stricted them only to the first half cycle, in order to have 
a closer look at quantum interference effects. Further- 
more, the uniform approximation employed in this work 
requires that we deal with pairs of trajectories (for de- 
tails see j an d such effects are expected to be related 
to single trajectories in a pair These issues were not 
the main objective of this work, and are open for further 
studies. 
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